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Synchronization effects in the dynamical behavior of elevators 
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We simulate the dynamical behavior of M elevators serving floors of a building in which a 
Poisson distribution of persons call elevators. Our simulation reproduces the jamming effect typically 
seen in large buildings when a large number of persons decide to leave the building simultaneously. 
The collective behavior of the elevators involves characteristics similar to those observed in systems 
of coupled oscillators. In addition, there is an apparently rule-free critical population density above 
which elevators start to arrive synchronously at the ground floor. 



The macroscopic effect of collective synchronization 
is a remarkable phenomenon occurring in several fields 
and is presently a topic of very active research. Effects 
of collective synchronization can be observed in situa- 
tions as diverse as, for example, the dynamics of sliding 
charge-density waves [l|, the phase-locking of relativis- 
tic magnetrons , and many other situations of interest 
in physics and engineering In addition, as early ob- 
served by Winfree [Q, synchronization is a phenomenon 
that also occurs in virtually all levels of biological organi- 
zation. Interesting examples of mutual synchronization 
involving biological phenomena are epileptic seizures in 
the brain, electric synchrony among cardiac pacemaker 
cells, flashing synchrony among swarms of fireflies, crick- 
ets chirping in unison, synchronization of menstrual cy- 
cles in groups of women, and several others ||^. A com- 
mon way of investigating the dynamics of such collections 
of oscillators is by assuming a relatively weak coupling 
K among them together with a random distribution of 
eigenfrequencies. Increasing coupling from zero one finds 
incoherent collective behavior up to a critical threshold 
after which the system starts to behave in a state dis- 
playing relatively high degrees of synchronization among 
the individual oscillators composing the system. A con- 
venient way of dealing analytically with several coupled 
oscillators is via a model proposed by Kuramoto Q , who 
assumed the dynamics to be ruled by the set of equations 
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i = l,2,...,iV. (1) 



In this equation N is the number of coupled oscillators 
while Oi and LOi are their phase and natural frequency, re- 
spectively. Recent progress in the understanding of the 
dynamical behavior of systems containing relatively few 
degrees of freedom along with a strong motivation to un- 
derstand nontrivial temporal and spatial structures has 
led to much interest in models for spatially extended sys- 
tems 0. In several of these investigations, not only the 
somewhat more traditional mathematical approach of us- 



ing sets of either ordinary or partial differential equations 
has been used, but new aspects have been uncovered by 
also using coupled map lattices and cellular automata 0] 

The purpose of this paper is to report synchronization 
effects observed in a computer simulation of the dynami- 
cal behavior of M elevators serving N floors of a building. 
Specifically, we report here the dynamics of the Feier- 
abend ejfect, that is, the jamming effect occurring at a 
certain time when a large number of occupants of a build- 
ing decide to leave it nearly simultaneously. Macroscopic 
synchronization effects between elevators were observed 
under a number of different conditions, all of them varia- 
tions around a simple set of "working rules" assumed to 
define the dynamical behavior of the elevators. These 
rules were chosen in such a way to provide a reason- 
ably fair representation of what one usually gets from 
not very sophisticated, not computerized elevators. The 
definition a set of rules controlling the behavior of eleva- 
tors is not a completely trivial task since the movement 
of elevators depends on a quite large number of factors. 
To start, while one might easily agree on a working def- 
inition for small and large buildings, it is hard to dif- 
ferentiate unambiguously the "transition" from small to 
large buildings. This question is relevant because eleva- 
tor rules good for small buildings are not practical for 
higher buildings and vice- versa. In higher buildings, it 
is common practice to divide the total number of eleva- 
tors into smaller groups in order to ensure good service 
for all floors. This division is not unique: one may use 
some elevators to serve only odd-numbered floors while 
the others serve even- numbered floors; or one may decide 
that a number of elevators will serve only the first, say, 
10 floors, the next will serve floors 10 to 20, and so on. 
It is not difficult to imagine a number of other possible 
combinations. In addition to all this, there are a number 
of specific details that need to be defined in order to op- 
timize the process, i.e. to minimize the time needed for 
a person to move between floors and to minimize energy 
consumption. 
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In what follows we will assume the ground floor to be 
the only possible way to exit the building and that all 
calls for elevators are from passengers wanting to exit 
by going to this floor. This choice reflects closely the 
phenomena that we want to consider. Each elevator can 
carry no more than 20 passengers simultaneously. The 
specific working rules of the elevators are assumed to be 
as follows. Whenever free elevators are available, every 
call originating from a floor from which there were no 
previous unresolved calls will set the closest free eleva- 
tor in motion. This will occur regardless of the relative 
position of the floor originating the call and of any pre- 
existing movement of elevators. After fulflUing the origi- 
nal call that set it moving, every elevator will always try 
to attend all other calls appearing in the floors lying be- 
low it, namely, it will stop at all levels lying below. While 
in practice one knows that an elevator already full can- 
not accept new passengers, real elevators are frequently 
not able to recognize and use the information that they 
are full to avoid stopping unnecessarily at lower floors. 
We simulate both situations: "naive" elevators, elevators 
which even after being full will continue to stop for ev- 
ery call from below, and "smart" elevators, elevators that 
once full will ignore all subsequent calls, moving straight 
to the lowest floor. Empty elevators arriving at empty 
floors, i.e. floors that have already been served, or on the 
ground floor will stop and either wait for a new call or 
move to floors where there are passengers calling, which 
could not be served earlier. 

The updating of the dynamics is done at regularly 
spaced time intervals, controlled by the discrete ticking of 
a clock. At every clock-step we assume p^{nk) new pas- 
sengers to call elevators from the k—th floor, distributed 
according to a Poisson law 

Pa.K) = A*"^^, (2) 

where /i is the Poisson parameter. At every clock-step 
the location and movement of every elevator will also be 
updated. The speed of the elevators is one floor per time- 
step, either up or down. Whenever an elevator stops in 
some floor, it remains in the floor during a few clock 
units. This time is intended to represent an average time 
consumed in opening the door, moving passengers and 
closing door. In the present simulations we assumed this 
time to be 5 clock units. In absence of any new call, 
elevators will preserve their previous state of movement. 
Some elevators may stop, if they arrive on the first floor. 
Otherwise, elevators will react to new calls as described 
above. To investigate the dynamical behavior as a func- 
tion of /i we perform simulations over relatively long time 
intervals (typically t = 5000 clock steps) accumulating 
several quantities at every clock-step, e.g. the total num- 
ber of elevators that arrive at the ground floor, the num- 
ber of times that L elevators arrive simultaneously at the 
ground level, the cumulative number of passengers trans- 
ported, the number of passengers waiting in each queue, 
etc. 



Based on these assumptions we wrote a FORTRAN 
program to simulate the dynamics of a few different 
buildings with several possible number of elevators. Sim- 
ulations were done over long time intervals and for differ- 
ent values of the Poisson parameter fi. Apart from stor- 
ing the aforementioned relevant quantities describing the 
process, the program could also be run interactively, gen- 
erating a PostScript on-screen animation of the dynam- 
ics, which could be easily video-recorded. An schematic 
snapshot of a video-frame is given in Fig. 
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FIG. 1. Schematic representation of a typical building. 
The black shading of each elevator is proportional to the num- 
ber of passengers that it carries. 

The shading of each individual elevator is proportional 
to the number of passengers that it carries. Further (not 
shown in the schematic figure) , the actual number of pas- 
sengers waiting on every floor was displayed on the left 
of the rectangle representing the building by a line seg- 
ment with length proportional to the population density. 
Line segments on the right side indicated the cumulative 
number of passengers transported since the beginning of 
every run. The majority of runs was done for M = 10 
and N — 50. For numerical values of this order, the po- 
sition, the population density inside every elevator (their 
shading) and the population density outside every floor 
(line segments) may be comfortably updated on screen at 
every time step, producing an actual movie of the dynam- 
ics. In this way we are always able to visualize simulta- 
neously the up-and-down movement of all elevators and 
to monitor the behavior of the population inside every 
elevator and on the several waiting queues. By monitor- 
ing the number of passengers in the waiting queues and 
inside elevators, it was possible to recognize many inter- 
esting effects such as, for example, when the number of 
elevators was well adapted to the size of a building, to the 
population density living in the building, etc. Simulation 
of a particular building consisted essentially of "playing 
the elevator game" for different values of /i and observ- 
ing the dynamics on screen while cumulating important 
numbers for later analysis. 
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Figure g shows for /i between and 0.1 the number of 
passengers transported and waiting elevators after 5000 
clock steps. One clearly sees the existence of an oper- 
ational "jamming threshold": for /i slightly above the 
critical value fic — 0.03 the 10 elevators are not able to 
cope with the traffic anymore: the number of waiting pas- 
sengers increases rapidly while the number of passengers 
transported tends rapidly to the maximum limit afforded 
by the capacity of the total number of elevators. 
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FIG. 2. Accumulated population density after 5000 time 
steps. Above /i ~ 0.03 the transport is not efficient anymore. 

Figure |^ shows a plot the critical values of /i obtained 
by simulation the dynamics of buildings containing 10 < 
TV < 80 floors and served by either M = 1, 2, 5, 10 or 20 
elevators. 




Number of floors 

FIG. 3. Critical values of /i as functions of the number 
of floors and number of elevators. As seen in Fig. 2, beyond 
Hcrit the transport is not efficient anymore. The characteristic 
exponent underlying all curves is ~ A*'"^'**. 



As one sees from the graph, the behavior on a log-log 
plot may be very well fitted by straight-lines. By chang- 
ing the size of the building, the capacity of the elevators 
and/or the number of elevators serving the building one 
obtains curves displaying qualitatively identical behav- 
iors. For example, for N — 50, we find ficrit — 0.00312-Af 
for 1 < M < 25. A typical recording of the time evo- 
lution of all elevators is shown in Fig. ^ which displays 
simultaneously the movement of all 10 elevators as a func- 
tion of time. Each line in this graph represents the ac- 
tual trajectory followed by an individual elevator. The 
10 lines are not individually discernible in Fig. ^: for 
/i = 0.01 because the 10 elevators move in an uncorre- 
lated way with their individual trajectories significantly 
overlapping each other, and for higher fi values because 
they then synchronize. 
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FIG. 4. Up-and-down movement of 10 elevators as a func- 
tion of time for three different Poisson distributions of calls. 



The point of the figure is to show that for /i— values 
of the order of ~ 0.1 the elevators move always stochas- 
tically, i.e. independently of each other, refiecting essen- 
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tially the randomness of the arrival of calhng passengers 
in each floor. As n increases one sees the appearance 
of white intervals along the lowest level. Such intervals 
represent quiescent periods, when all elevators are busy 
elsewhere. Thus, as /i increases there are sequences of 
roughly periodic time intervals where the majority of el- 
evators bunch at the lowest quota (ground level), indi- 
cating a relatively high synchronization of the arrivals 
on the ground floor as the number of calls increases. 

An alternative way to recognize the synchronization is 
by monitoring what happens during the time evolution 
of the cumulative number of elevators arriving on the 
ground floor. This is recorded in Fig. ^. 




FIG. 5. Time evolution of the cumulative number of ele- 
vators arriving on the ground floor. The lower picture is a 
magniflcation of the portion inside the rectangle in the upper 
one. 



events recorded. As fi increases, there is a redistribution, 
with many more elevators arriving simultaneously. For 
/i = 0.03, near the jamming threshold, one may already 
recognize visually the presence of an underlying period- 
icity in the number of simultaneous arrivals. 
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Histograms showing the number of simultaneous 
a function of time and population density. 



For small population densities one observes essentially 
a linear increase. For higher densities the number of ar- 
rivals increases discontinuously, roughly periodically, via 
"quantum" jumps. These increases occur in the rela- 
tively short time intervals when several elevators arrive 
simultaneously. The smallest (and barely visible) vertical 
jump corresponding to the arrival of a single elevator may 
be more easily recognized from the curve corresponding 
to the lowest density, i.e. /i = 0.01. The synchroniza- 
tion and its corresponding characteristic frequency may 
be also recognized from Fig. ^. As a function of time, 
this figure shows the number of coincidental arrivals at 
the ground floor for three different values of /i. For low /i, 
the most salient feature is that elevators tend to arrived 
independently, as indicated by the large number of single 



With respect to Fig. ^ it is important to notice that the 
figure gives the total number of elevators arriving at ev- 
ery single clock step. This is a very stringent definition of 
synchronization because it only considers as "simultane- 
ous" those arrivals happening at exactly one single time 
step. A more tolerant characterization would very likely 
involve considering as "simultaneous" all arrivals occur- 
ring within a time- window larger than a single clock unit. 
The increase in the number of quasi-simultaneous arrivals 
seen in Fig. ^ is responsible for the regularly spaced jumps 
appearing in Fig. 

To recognize quantitatively the frequencies underlying 
the regularities seen in the spacings of Fig. ^ we com- 
puted Fourier transforms for the time series shown in 
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Fig. |5[ Three such transforms are given in Fig. 0, for 
values of /i as indicated. 
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FIG. 7. Fourier transform for three of the time series given 
in Fig. 5, as indicated by the value of fi. 

From these figures one recognizes the appearance of a 
clear peak as fi increases, giving the characteristic fre- 
quency with which elevators synchronize among them- 
selves. 

So far, all results presented were obtained assuming 
"naive" elevators, i.e. elevators that while going down 
will stop for every call originating from floors below them, 
even after being full. Although this is a quite common 
behavior in real elevators, it is of interest to investigate 
what happens if elevators are smart enough to notice 
when they are full and use this information to avoid un- 
necessary stops. This question is important because full 
elevators will not further contribute to the transport and 
their excessive stops will only lead to a decrease of the 
average speed and to a bunching at the highest elevator 
velocity (versus floor coordinate) gradient. ^From these 
heuristic considerations one may expect synchronization 



to disappear for smart non-stopping elevators. Figure g 
shows the evolution of the total number of passengers 
waiting for elevators as a function of /i for both types of 
elevators. 
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FIG. 8. Total number of waiting passengers for both naive 
and smart elevators. Numbers refer to the elapsed time after 
which measurements were done, in units of 1000 clock-steps. 

The four sets of curves display the total of passen- 
gers waiting for elevators after the indicated number of 
clock steps, divided by 1000. It is clear that beyond 
fj^crit the transport is not efficient anymore and the num- 
ber of passengers will diverge if one waits long enough. 
These curves are given to indicated the how this diver- 
gence sets in four both types of elevators. Notice the 
differences in the behavior around fj-crit- The differences 
may be also recognized from Fig. ^, which shows on a 
single plot two sets of curves: the total number of pas- 
sengers transported after 5000 clock units and the total 
number of arrivals at the bottom of the building. 
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FIG. 9. Total number of transported passengers and total 
number of arrivals at the ground floor for both rules, naive 
and smart, measured after 5000 clock steps. 

As expected, one sees that smart elevators have larger 
^crit, being therefore able to provide a much more efB- 
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cient transport. We also made plots similar to Fig. ^ to 
check that synchronization is still happening. The figure 
obtained is similar to Fig. ^ with smart elevators show 
synchronization over a very small fj, interval. 

One interesting aspect of the simulation is the way 
in which the jamming of passengers occurs. While the 
asymptotic number of transported passengers is clearly 
limited by the number of elevators and their capacity, it 
is somewhat surprising to find that as /i increases, this 
limit is attained after an "overshoot" of the number of 
transported passengers occurring exactly at the jamming 
threshold. As clearly seen from Fig. ||, rather than con- 
verging monotonically to the saturation level dictated by 
the number and capacity of the elevators, the most ef- 
ficient transport occurs around the jamming threshold 
where one finds a marked discontinuity in the derivative 
of the number of transported passengers as a function of 
/i. Beyond the jamming threshold the simulations indi- 
cate the existence of two basic behaviors while passengers 
start to accumulate in all floors, waiting for elevators that 
tend more and more to appear full. Immediately after the 
jamming threshold there is a "relaxation interval" of fi 
values characterized by negative derivatives. After this 
regime the number of transported passengers is indepen- 
dent of fi. All in all, the best compromise from the point 
of view of the passengers wanting to leave the building is 
to plan the building to operate at the jamming threshold. 
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